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Abstract: We investigate the primordial non-Gaussianities from the trispectra in multi-field inflation models, which can 
be seen as generalization of multi-field fc-inflation and multi-DBI inflation. We derive the full fourth-order perturbation 
action for the inflaton fields and evaluate the four-point correlation functions for the perturbations in the limit ^ 1 
and Ce <C 1. There are three types of momentum-dependent shape functions which arise from three types of four-point 
interaction vertices. The final trispectrum of the curvature perturbation can be expressed in terms of the deformations and 
permutations of these three shape functions, and is determined by Ca, Ce, A , 11 which depend on the non-linear structure of 
the model and also the transfer function Tus- We also discuss the parameter space for the trispectrum and plot the shape 
diagrams for the trispectrum both for visualization and for distinguishing different shapes from each other. 
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^ Explicit Expressions for {QaQaQsQs)i 



1. Introduction 

Current observational data support cosmological inflation greatly in which primordial perturbations assumed respon- 
sible for Cosmic Microwave Background anisotropics and Large-Scale Structure formation are generated from quantum 
fluctuations and stretched to superhorizon scales during inflation (see e.g. ^ for a review). Actually, inflation itself 
is not a single model, but rather a theoretical framework. One of the most robust predictions of inflation is the almost 
scale-invariant, Gaussian and adiabatic primordial fluctuation. However, from the point of view of power spectrum (i.e., 
Fourier transformation of the tree-level two-point equal-time correlation function) of the cosmological perturbation, many 
inflation models are "degenerate". The theory and observation of power spectrum do not give us an unique theory of 
inflation. Phenomena beyond linear-order have been extensively investigated over the past several years. 

The most significant progress beyond power spectrum is the investigation of statistical non-Gaussianities of CMB 
anisotropics and primordial fluctuations (see e.g. for a nice review). From the field theoretical point of view, 
non-Gaussianity describes interactions among perturbations, which will cause non-vanishing higher-order correlation 
functions. Such interactions are mandatory in any realistic inflationary models, which come from both the non-linear 
nature of gravitation and the self-interactions of inflation field(s). In standard slow-roll inflation scenario, however, the 
non-Gaussianities have been proved too small to be detected (see also [P), even with PLANCK satellite [10|. Thus, 



any detection of non-Gaussianities would not only rule out the simplest inflation models but also give us valuable insight 
into fundamental physics of inflation l|l||l||l|,|l|,|l||l|,|l7[|l8[|l||. 

Generally speaking, primordial non-Gaussianities can be large, if one or more of the following four conditions are 
violated: 1) slow -roll conditions, 2) canonical kinetic terms, 3) single field and 4) Bunch-Davies vacuum [^. For single- 
field models, for example, various possibilities have been investigated in order to generate large non-Gaussianities by 
introducing complicated kinetic terms pO| , pT| , p2[ p4| , p5[ ^ , 27 28], ^ ^,31 3^ which belong to the more general 



class of fc-inflation models in which the small-speed of sound Cs ^ 1 enhances the derivative-coupling of 

perturbations, or other mechanisms which enhance the interactions during inflation or non-linearities during superhorizon 
evolution [|| g IH I?! g ||, |^ 0, H . 

From the point of view of data analysis, there are two typical types of non-Gaussianities which are most interesting: 
the so-called "local" type and the "equilateral" type. The former describes the strength of modulations of short wavelength 
perturbation modes by long wavelength modes, while the later describes the correlations among modes with similar wave- 
lengthes. It has been clear from the studies of non-Gaussianities in single-field models that, single-field inflationary models 
with non-canonical kinetic terms can have large "equilateral"-type non-Gaussianities, since the "derivative-coupled" inter- 
actions among different modes are enhanced by the non-canonical structure of the kinetic term. While single-field models 
can never generate large "locaf'-type non-Gaussianities, unless the generalized slow-roll conditions are abandoned. 

Multiple field inflationary models provide us with more possibilities to generate large non-Gaussianities during in- 
flation [|7[|8|,|9l|^|l|||||||,|5|,||g|7l|||9l00||0||,|6[|3]. It has been shown fliat multi- 



inflaton models with non-canonical kinetic terms are also mainly characterized by equilateral-type non-Gaussianities 
rather than local-type. Local-type non-Gaussianities can be generated during superhorizon evolution of cosmologi- 
cal perturbations, for example in inhomogeneous "end-of-inflation" models [ |75| , [T^ fl\ [79[ [78[ | or curvaton scenario 
[p4ip5ip6lp7i^p9l |l00l|l0lt . 



In this work, we investigate primordial non-Gaussianities from the trispectra (i.e. Fourier transformation of equal- 
time four-point correlation functions) of cosmological perturbations in general multiple field models. We consider the 
model with general scalar-fields Lagrangian of the form P{X" , cj)^). This form of La grangian was first investigated in 

p2| . It includes single-field models with non-canonical kinetic term [ p4| , , multi-field k-inflation [ |73| , pO| and multi- 
field DBI models [||||,|lH|] as special cases, and thus deserves detailed studies. Bispectrum in general multi-field 
models with small speed of sound has been investigated in ^% ^ |60[ p^ ]. Trispectra in single-field models have been 
investigated by various author [|2[ |6[|^ ||, |9[ |3|, |l], |2[ and in multi-field models until very recently |4[ ||, m . 



We start from the the full fourth-order action for the scalar perturbations and then evaluate the leading-order four-point 
functions around sound-horizon crossing in the approximation Ca ^ 1 and Ce ^ 1, where Ca and Cg are propagation speeds 
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of adiabatic and entropy modes respectively. In this work we restrict the calculations to the contributions from the so-called 
"contact interaction", i.e. the contributions from the "four-point direct interaction vertices"'. As in the case in general 
single-field models, in the leading-order, we find three fundamental shape functions Ai(ki, fc2, ^3, ^4), A2{ki, fc2; fci, ^2) 
and A3{ki, k2\ fcs, k^), corresponding to three different types of four-point interaction vertices. Moreover, in addition 
to the pure adiabatic four-point function {QcrQaQaQa)^ there also exist pure entropy four-point function^ {QsQsQsQs) 
and one mixed contribution {QaQaQsQs) (we denote Qo- and Qs as the adiabatic and entropy modes respectively). The 
four-point correlation functions (QaQaQaQa), {QaQaQsQs) and (QsQsQsQs) can be generally written in terms of 
various deformations and permutations of these three types of shape functions. Since the observable is the curvature 
perturbation TZ, we also investigate the trispectrum of the curvature perturbation (TZTZTZTZ) on large scales. 

In this work, we assume the propagation speeds of adiabatic mode and entropy modes be different: Ca 7^ Ce. Thus, 
unlike the multi-field DBI models |8[ ^ where Ca — Ce = Cs and we can abstract common shape functions for 
both (QaQaQa) and (QaQsQs), in this work we will see that Ca, Ce enter the definition of the shape function intrinsically. 
More precisely, as we will see, the trispectrum of the curvature perturbation is determined by four parameters Ca, Ce, A, 11 
which arise from the non-linear structure of the theory and also T-jig which characterizes the transfer from entropy modes 
to adiabatic modes during superhorizon evolution. 

This paper is organized as follows. In section 2, we describe the model and set up the perturbation theory in spatially- 
flat gauge, then we briefly review the linear perturbation and the power spectra in section 3. In section 4, we give the 
full fourth-order perturbation action and derive the interaction Hamiltonian in the interaction picture. And we calculate 
the four-point functions and the trispectra in section 5. In section 6, we give a general discussion on characterizing the 
trispectrum and we plot the shape diagrams and investigate the non-linear parameters. In the end we make a conclusion 
and discuss some related issues. 

In this paper, we choose c = h = Mpi = = 1. 



2. Basic Setup 

2.1 Model and Background 

In this work we consider a general class of multi-field models containing scalar fields coupled to Einstein gravity. The 
action takes the form 

5 = 1 d^xV^\^ + P{x",^')\ , (2.1) 

where (/ = 1, 2, • • • , A/^) are scalar fields acting as inflaton fields, and 

X^'^ = (2.2) 

is the kinetic term (matrix), is the spacetime metric tensor with signature (— , +, +, +). "/, J"-indices are raised, 
lowered and contracted by A/^-dimensional field-space metric Gjj = Gjj{(j)'). This form of Lagrangian includes multi- 
field k- inflation and multi-DBI models as special cases^. 

Modern "action approach" of cosmological perturbations is based on the ADM formalism of gravitation, in which 
spacetime metric is written as 

ds^ = -N^dt^ + h,j{dx' + N'dt){dx^ + N^dt) , (2.3) 

with N — N{t, x) is the lapse function and Ni — Ni{t, x) is the shift vector, hij is the spatial metric on constant time 
hypersurfaces. The ADM formalism is convenient because the equations of motion for N and are exactly the energy 



and momentum constraints which are easy to solve. Under the ADM formalism, the action (2.1) can be written as (up to 
total derivative terms) 

J dtd^xVhN Q^'^^ + ^ {EtjE'^ -E^)^ + J dtd^xVhNP, (2.4) 



'Obviously the tree-level four-point coiTelation functions have two origins: one is the four-point vertices, the other is con'elating two thi'ee-point 
vertices, e.g. with scalar modes 0,^] or graviton [po|. A full analysis should include these two contributions together. 

^This is different from the case of bispectrum. There is no leading-order pure entropy three-point correlation (QsQsQs) and the mixed contribution 
is of the form (QaQsQs), see |^ ^ for details. 

^ For example, multi-field k-inflation has the scalar-field Lagrangian as P{X,(j>^), where X = tiX^'^ = GjjX^'^. While in multi-field DBI 
models, P{X''',4,') = -jfjT^ (yv- l) - V{4i') with© = 1 - 2fGijX'-' + Apx''jXf - 8f'-^x[' XjX^'^ + le/^xj^X^Xf X^'. 
See later discussion for details. 
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where h = det hij and the symmetric tensor 

£y = i - V,Nj - VjTV,) , (2.5) 

with Vi is the spatial covariant derivative defined with spatial metric hij and E = tvEij = h^^ Eij. R^^^^ is the three- 
dimensional Ricci scalar which is computed from the spatial metric hij . In ADM formalism, spatial indices are raised 
and lowered using hij and /i*^ . 

In ADM formalism, the kinetic matrix X^'^ can be written as 

X'' = -\h^'d^<l^'d,^' + i^v'v' , (2.6) 

where 

= (p' - N'VicI)^ . (2.7) 

2.1.1 Equations of Motion 

The equations of motion for the scalar fields are 

V4P(/,/)9'^0O+^../ = O' (2-8) 
where is the four-dimensional covariant derivative. Here and in what follows, we denote 

_ dP _ d^P 

P,iiJ) = ' P,(i.J)(KL) = -qxTTqxkl ' (2.9) 

for shorthand. 

The equations of motion for N and Ni are Hamiltonian and momentum constraints respectively. 



V 



(2.10) 



2.1.2 Background 



In this work, we investigate scalar perturbations around a flat FRW background, the background spacetime metric takes 
the form 

ds2 ^ _^^2 ^ a'^{t)5,jdx'dx' , (2.11) 
where a{t) is the so-called scale-factor The Friedmann equation and the continuity equation are 

3 3 ^ '^-^''^ ' ' (2.12) 

p=-iH{p + P). 

In the above equations, all quantities are background values. From the above two equations we can also get another 
convenient equation 

TT vIJ n 

The background equations of motion for the scalar fields are 



H = ~X''P,uj). (2.13) 



P{iJ)4>' + (iHP,(ij^ + P,(ij)) - P,j = , (2.14) 



where Pj denotes derivative of P with respect to cj) : Pj 

In this work, we investigate cosmological perturbations during an exponential inflation period. Thus, from ( ^.13[ ) it 
is convenient to define a slow-roll parameter 

e^-A.^MM (2 15) 
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2.2 Perturbation Theory in Spatially-flat Gauge 

The scalar metric fluctuations about our background can be written as (see ||] for nice review of the theory of cosmo- 
logical perturbations) 



SN = a, 



(2.16) 



where a, /?, i/" and E are functions of space and time'*. The scalar field perturbations are denoted by Scf)^ = . 

Before proceeding, we would like to analyze the (scalar) dynamical degrees of freedom in our system. In the begin- 
ning we have A/^ + 4 apparent scalar degrees of freedom. The diffeomorphism of Einstein gravity eliminates two of them^, 
leaving us A/^ + 2 scalar degrees of freedom. Furthermore, two of these M + 2 degrees of freedom are non-dynamical. 
In ADM formalism, these are just the fluctuations 5N = a and 6Ni — dif3. Thus, there are J\f propagating degrees of 
freedom in our system. As has been addressed, the diffeomorphism invariance allows us to choose convenient gauges 
to eliminate two degrees of freedom. In single-field models, there are two convenient gauge choices: comoving gauge 
corresponding to choosing Sifi = E = or spatially-flat gauge corresponding to ip — E ^ 0. In multi-field case, the 
comoving gauge loses its convenience since we cannot set (5/9 = in multi-field case. Thus, in this work we use the 
spatially-flat gauge^. 

In the spatially-flat gauge, propagating degrees of freedom for scalar perturbations are the inflaton field perturba- 
tions Q^{t, x), while SN and SNi are non-dynamical constraints. In this work, we focus on scalar perturbations^. The 
perturbations take the form 



K 

1 



N 



(2.17) 



ai + a2 



72i 



where (/)o(i) is the background value, and a„, Oni are of order 0{Q^). 

The next step is to solve the constraints Q!„, /3„ and 6'„j in terms of Q^. Fortunately, in order to expand the action to 
fourth-order in , the solutions for the constraints up to the second-order are adequate. In general, we need the solutions 
for constraints N and Ni up to ©(QI"/^!) order ([n/2] means the largest integer < n/2), if we want to expand the action 
to 0{Q'^) order and to calculate the n-point correlation functions. 

2.2.1 Solving the Constraints 



At the first-order in , a particular solution for equations (2.10) is: 

1 „ 



ai 



2H 



?H = . 



(IJ) 



2X'^^P 



{IJ){KL) 



X 



IJ 



H 



-P 



(KL) 



■Q 



(2.18) 



-3HP^^ij^<p'Q' ~ P^dj^KQ^'^X'-' + PjQ' 



Here and in what follows, repeated lower indices are contracted using 5ij, and 
should be understood in fourier space. 



didi. 9 is a formal notation and 



''This form of ansatz corresponds to 5goo = 1 ~ + NiN' and Sgoi = Ni. 
^See Ipti for a detailed discussion on the gauge issue of cosmological perturbations. 



*In the spatially-flat gauge, it is obvious that the dynamical degrees of freedom come from the scalar fields peiturbations S4>^ . In Einstein gravity, 
the seal ar-tvpe metric perturbations are essentially non-dynamical. While in non-Einstein gravity, e.g. in recently proposed Hofava gravity, this may not 
be the case 

^In general, it is well-known that in the higher-order perturbation theories, scalar/vector/tensor perturbation modes are coupled together However, 
from the point of view of the perturbation action approach, t hese couplings are equivalent to exchanging various modes. For example, the effects of 
tensor perturbation on scalar modes have been investigated in [ 11 3 1 through graviton exchange approach. In this work, we focus on interactions of scalar 
modes themselves, thus we can neglect tensor perturbations. 
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Similarly, at the second-order, we have 

ZH 



with 



djai 



P 



,{IJ){KL) 



and 



n = 



P.{IJ) - P.{IJ){KL)<Pa <Po 



2a2 



+ {P,{KL)M - ^'oHP{IJ){KL)m) Q^' (^f Q''^ - «1 'Z'^^^) 

+ \ [P.KL - H^P^iLDKL) Q^'Q'' + 2ai [P,{IJ){KL) (^f Q^^' " «! ^^'^o) + ^,UJ>ifQ 

- P,{IJ){KL) {'pi'^Q''^ - «1 ^Uo) + P,(IJ)kQ'' - 2aiP(/J)o] '^i'i'Q'^ 
-P{iJ) (q'Q' -2^\^Nld,Q'^ 



+ 2 [P,{KL)(MN) - 'i'o<PoP,{IJ){KL){MN) 



(2.19) 



(2.20) 



(2.21) 

It is useful to compare our above results with previously known results ^ ^ ^], which the above results reduce 

to in various limits. 



3. Linear Perturbations 



3.1 Instantaneous Adiabatic/entropy Modes Decomposition 

In multi-field inflation models, it is convenient to decompose perturbations into instantaneous adiabatic and entropy 
perturbations [ ]80| , pl| ]. The "adiabatic direction" corresponds to the direction of the "background inflaton velocity" 



ei 



where we define a 



PijK)<i>'<i>^ " 



o- = yP,{jK)(l}-^(l> 



(3.1) 



(3.2) 



which is the generalization of the background inflaton velocity. Actually & is essentially a short notation and has nothing 
to do with any concrete field. Note that a is related to the slow-roll parameter e as ct^ — 2H^e. 

We introduce (JV — 1) basis e^, (n — 2, • • • ,J\f) which are orthogonal with e{ and also with each other. The 
orthogonal condition can be defined as 

P,{IJ)^ii^ri = ^mn ■ (3.3) 

Thus the scalar-field perturbation can be decomposed into instantaneous adiabatic/entropy basis; 

Q^ = e^Q™, m^l,---M. (3.4) 
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Up to now our discussion is rather general, without further restriction on the structure of P{X''^ , cj)^). In this work, 
we consider a general class of two-field model, with the Lagrangian for the scalar fields as the following form^: 

P{X",<j>') = P{X,Y,(f>'), (3.5) 

with X = X\ ~ GijX^'' and Y = XjX/. This form of Lagrangian not only is the most general Lagrangian for 
two-field models and thus deserves detailed investigations, but also can make our discussions on the non-Gaussianities in 
two-field models in a more general background. For two-field case, higher contractions among Xjj can be expressed in 
terms of X and Y, e.g., 

X'jX-^X^ ^-^ + ^XY, 

x'jXix^xf ^-\x^ + x^Y + ir^ ^ 



etc. The model ( 3.5 ) includes multi-field fc-inflation and two-field DBI model as special cases. For example, in multi-DBI 
model the Lagrangian is P = — j^^j-) f — l~j — V"(<^^) with 



V = det (gS - 2/x5) 



(3.6) 



= 1 - 2fGijX'-' + Afx'fxf - S/^xf X/X^'l + IQf^X^lxjX^X^ 

This expression for determinant T) is general. In this work, we focus on two-field case, thus the last two terms exactly 
vanish, leaving as effectively P = l-2fGijX" + 'ipx\'xj\ In terms of this is P = l-2fX + 2p {X^ - Y). 
At this point, it is convenient to introduce two parameters: 

A = X'^P.xx + "ix^P.xxx + 21 YPy + &Y^P.yy + Iy^P.yyy 



+ 4 {X^YP.xxY + 2XYP,xY + 2XY^P^ XYY) , 
n = X'^P^xxx + \x^P,xxxx + ^ {21Y^P.YY + MY^P^YYY + ^Y^P.yyyy) (3.7) 

+ -rrX^Y Pxxxy + -^X^YP xxy + -rrX'^Y'^P xxyy + ^XYP_xy 
5 5 5 

H — ^XY'^Pxyy + -;rXY^P xyyy , 
5 5 



where all quantities are background values, and we have used Y — X^. As we will see later, although the X, In- 
dependences of P{X, Y, (f>^) in general can be complicated, the non-linear structures of P affect the trispectra through the 
above specific combinations of derivatives of P. 

3.2 Quantization and the Power Spectra 

After instantaneous adiabatic/entropy modes decomposition, at the leading-order, the second-order action for the pertur- 
bations takes the form' 



S2 ' — / dtd Xa ( ^JCmnQmQn ^^Smn9iQmdiQn ) ; (3.8) 



2 2a^ 
with 



(3.9) 

^ - 1 ) SlmSin + ( ^ - 1 ) {Smn - SimSin) , 



^This form of Lagrangian if motivated from that, for multi-field fc-inflation models |^^^], the Lagrangian is simply P{X, 4>^), In [ p2| | a special 

form of Lagrangian p(y, (t>') with Y = X + (x2 - XijX'-') was chosen in the in vest igation of bispectrua in two-field models, which is 
motivated by multi-field DBI action. In this work, we use the more general form of Lagrangian ( |3.5| ). 

'in ( [3.8| ) we neglect the mass-square terms as MmnQ^Q^ and the friction terms such as ~ QmQn- In general these terms may become important, 
especially they may cause non-vanishing cross-con'elations between adiabatic mode and entropy mode around horizon-crossing. See e.g. [^, [73[ 
for details. 
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where we introduce 



10 



2 Px + 2XPy 

c, = 



P.x + 2X [Pxx + 'iXPxY + 3P,Y + ^X^P.yy) ' 

p (3.1U) 



which are the propagation speeds of adiabatic perturbation and entropy perturbation respectively. It is useful to note that 
K-mn is diagonal, K-n — /C22 = 1/Cg and K.12 = ^21 = 0, as a sequence of adiabatic/entropy decomposition. 

Ca ^ Ce is a generic feature in multi-field models, this can also be seen explicitly from the definitions in ( 3.10 ), the speed 
of sound for adiabatic mode and entropy mode(s) have different dependence of P-derivatives". 

Obviously, and Qs themselves are not properly normahzed for canonical quantization. We may introduce new 
variables 

Qa=-Qa. Qs = -Qs, (3.11) 

which are canonically normalized variables, since after changing into comoving time defined by dt = adrj, the quadratic 
action takes the form 



52 = 



(3.12) 



j drid\ i \q'^ + {re + W) Ql - clidQ^f + + {re + W) Ql - clidQsf ■ 

The action ( |3.8| ) or ( 3.12 ) describes a free theory, which is easy to quantize. In canonical quantization, we write 

Qa{k,r]) = akUkiv) + a'Li^uKr]) , Qs{k,v) = akVkiv) + a^-k'^liv) , (3.13) 

where Uk {rj) and Vk iv) mode functions, which satisfy the corresponding classical equations of motion 

< + [cafc' - + ^')] uk^o, i'l + {cik^ - {re + m)\ = , (3.14) 

and can be easily solved in de Sitter approximation (a(ry) — — -fj^): 

Note that the mode functions are chosen so that when the modes are deep in the sound horizon, or equivalently 77 —00, 
they behave as free harmonic oscillators in Minkowski spacetime, i.e. 

nk{n) 71:1 ' ^'^"^^ vif • 

Moreover, Uk {rj) and Vk (rj) are normalized with Wronskian 

'Uk{v)uk{v) ~ Kiv)u'k{v) = Vk{'n)vk{v) - vl{ri)v'i.{ri) = i, (3.16) 

which is the condition for canonical quantization. 

Finally, what we are interested in are the tree-level two-point functions for Qa and Qs, defined as 

{Qa{kum)QAk2,m)) = {27:)^s\ki + k2)Gk, {m.m) , 
{Qs{ki,iu)Qs{k2,m)) = (2^)^(52(^1 + fc2)Pfei {m.m) , 

with 

Gk{vi,V2) = Uk{m)uUV2), Fk{vi,V2) = Vk{m)K{V2), (3.18) 

where Uk {rj) and Vk {f]) are the mode functions for adiabatic perturbation and entropy perturbation respectively: 



(3.17) 



(3.19) 



'"We use Ca and Ce rather than Co- and Cs in order to avoid possible confusion, since in the literatures Cs has special meaning, i.e. the speed of sound 
of perturbation in single-field models. 

"This fact was first point out apparently in js^ , p3| in the investigation of brane inflation models. See also ^ |8^, ^ [s^, ^(^ for extensive 
investigations on general multi-field models with different Ca and Cc. 
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The so-called "power spectra" for adiabatic perturbation and entropy perturbation are defined as Pa{k) = Gkirj^., rj^,) 
and Ps{k) = Fk{ri^,ri^), where 77, can be chosen as the time when the modes cross the sound-horizon, i.e. at c^k = aH 
for adiabatic mode and c^k = aH for entropy mode(s)'^. We have 



Pa*{k) 



2c,k^ ' 



Ps* (/c) 



(3.20) 



In the so-called comoving gauge, the perturbation is directly related to the three-dimensional curvature of the 
constant time space-like slices. This gives the gauge-invariant quantity referred to as the "comoving curvature perturba- 
tion": 



H 



Qa, 



(3.21) 



where a is defined in (3.2 1. The entropy perturbation Qs is automatically gauge-invariant by construction. It is also 



convenient to introduce a renormahzed "isocurvature perturbation" defined as 

S=-Qs- 
a 



(3.22) 



In the cosmological context, it is also convenient to define the dimensionless power spectra for comoving curvature 
perturbation and isocurvature perturbation respectively: 



i72 ^2 ^3 1 



IP 
a2 



2eca V27r 



ct2 27r2^ 
H'^ k^ 

Vs. = ^—Vs. = —-^Ps.{k) = — ( ^ 



1 

2ece 



(3.23) 



In the above results, all quantities are evaluated around the sound-horizon crossing. V-jz in ( 3.23 ) recovers the well-known 
result for single-field models 21 1 . In the case when Ca = Ce, the above results reduce to those in multi-field DBI 

model which has been investigated in |8 62 1. 



3.3 Superhorizon Evolution 

Actually, the inflaton fields perturbation, or Qc and Qs themselves are not directly observable. What we are interested 
in is the curvature perturbation. In single-field inflation models, this has no particular difficulties since the comoving 



curvature perturbation TZ is conserved on superhorizon scales [102, 103, 104|. Thus, it is sufficient to evaluate the 



correlation functions for i c- the curvature perturbation on uniform density hypersurfaces which coincides with —TZ on 
superhorizon scales, at the time of horizon-crossing. 

However, in contrast with the single-field models, in multi-field models, the curvature perturbation in general evolves 



after the horizon-crossing [107| (see also [109|). This is due to the fact that, in superhorizon scales adiabatic perturbation 



can be sourced by entropy perturbation(s) and there is a transfer between adiabatic/entropy modes. This can be clearly 
seen if we take the time derivative of the curvature perturbation [|80|] (see [O] for a recent investigation), in our model it is 



H a-^ 



where 



~ (l + c2)P. 



(3.24) 



(3.25) 



Due to the presence of isocurvature perturbation S, even on superhorizon scales \c^k/a\ <C 1, Tt w ^5 7^ 0, that is TZ will 
evolve. In general, on superhorizon scales, the evolution of curvature/isocurvature perturbations can be approximately 
described by 

TZ « aHS , S « (3HS . (3.26) 



'^In general multi-field models, adiabatic mode and entropy modes(s) with the same comoving wavenumber k exit the sound-horizon at different 
time, due to their different speeds of sound, ^ c^. This phenomena may cause subtle problems, such as the problem of decoherences. In this work, 
we neglect these problems. 
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The above equations have formal solutions 



with 



!ldt' P(t')H(t') 



1 Tns\ 
Tss) UA ' 



dt'a{t')Tss{t',U)H{t'), 



(3.27) 



(3.28) 



where is the time of horizon-crossing and t is some later time. Thus on superhorizon scales, the (time-dependent) 
power spectra for the curvature perturbation, isocurvature perturbation and also the cross-correlation between the two can 
be formally expressed as 

'Ps{t)^Tls{t,t.)Vs*, (3.29) 
Cns{t) = (nS) - TT^sit,U)Tss{t,U)Vs. ■ 

4. Non-linear Perturbations 

In this section, we derive the fourth-order action for the perturbations and the fourth-order interaction Hamiltonian (in the 
interaction picture). In the next section, we evaluate the four-point correlation functions for the perturbations. 

4.1 Fourth-order Perturbation Action 



From ( |2.4| ), Taylor expansion gives the fourth-order perturbation action from the gravity sector 

S'i^ J dtcPx y |-6i/2 (^4 _ iala2 + al) + (-a? + 2aia2) (^^'a) 



AH 



1 r 



d,djPid,d,Pi-{d^(3iy 



(4.1) 



2ai 



J [{d^djf3i) {d^dj(32 + a(,02,)) - 52/31^2/32] + ^ [2 9.9^/32 ft^2j + d^92jd(,92j) 
The scalar field sector is more complicated, 

St = j dtd^x {Pi + aiPy. + ^2^2) , 



(4.2) 



where P„'s are the corresponding parts in the expansion of P which are order 0{Q'^) respectively, which we do not wrtie 
explicitly here for clarity, and can be found in Appendix ^ 



4.1.1 Fourtli-order Action at the Leading-order 

In this work, we focus on the dominant contributions to the non-Gaussianities from the trispectra for Ca <C 1 and Cg ^ 1. 
In the case of bispectrum, in this limit, contributions from the gravity sector 5*1 and the metric perturbations themselves 
ai, (3i can be neglected in the leading order H owever, this is no longer the case for the fourth-order calculation. The 
orders of various quantities can be read from ( 2.18 )-( 2.19 ) and are summarized in Tab.|l]. 





Pi 




(X2 


P2 


02i 








0{\lcl) 


Oil/cl) 





Table 1: Summary of the order of constraints. Cs denotes both Ca and Cc. 
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Moreover, in the slow -roll limit, the fourth-order action from the gravity sector can be approximately neglected, since 



— I ~ e. Thus, at leading-order in slow-roll, the fourth-order perturbation action reads 



Q(main) 



dt(fx 



' T3 \riJ\rl<^l^,'Tj \rlj ^KL \r M N 



1 ^ p yU vKL yMN yPQ 

+ -^K{IJ){KL){MN){PQ)-^1 ^1 ^1 



+P. 



1 



\{ij){KL)-^d,Q'm'd,Q''dj( 



(4.3) 



where 



^IJKL = -^P,{I.J){KL) + -^P,{IJ){MK){NL)4>0^ <i>0 + -^P,{MI){NJ){PK){QL)4'o' 'i'O 4'o 'i^o ' 



QlJKL = P,(IJ){KL) + P,{KL){MI){NJ)WPo , 

for simplicity. The explicit expression for X^'^ can be found in Appendix ^ 

After adiabatic/entropy decomposition, the fourth-order action (4.3) can be written as 



(4.4) 



(main) 



where 



dtcfix 



- mnpq 



mnpq QmQnQpQq ^ 2 ^fnnpq QmQn^iQp^iQq ^~ g 4 ^mripq ^iQrn^iQri^jQp^jQq 



\P,{IJ){KL) ei^eie^e^ + j (p^(pq)^^4>^^ P,{ij){MK){NL)ei' ei,eie^e 



MNP, 



e 

nr. 



mnpq 



^mnpq 

4.2 Interacting Hamiltonian 



P.(iJ){KL) ei^e^epBg 
P..(iJ){KL) eineiepBg 



P 



(4.5) 



(4.6) 



In the operator formalism of quantization, interaction Hamiltonian is needed. The interactions of cosmological pertur- 
bation are in general contain time derivatives, which are different form ordinary field theory where interactions are local 
products of fields. Thus, in order to get the corresponding Hamiltonian, we should use its definition 



(4.7) 



where C is the Lagrangian containing 2nd, 3rd and 4th-order terms. The 2nd-order part is given in (3.8). The 3rd-order 
part has been derived in 



SI 



'(main) 



dtd X 0, I ^^mnlQmQnQl 



1 



(4.8) 



with 



' — 'mnl 



PiMN)C^''4>'' 



P 



,{MN)\ 



IM IN 



P 



.{iK){jL){PQ)eie^e(e^^e^ei 



(4.9) 



From (3.8 ), (4-. 5 ) and (fl-.8|), through a straightforward but rather tedious calculation, we get 



7^4 — QmQnQpQq 
1 



)'^rs ^rmn ^spq 



- mnpq 



1 . 



QniQn^iQp^iQq \ 



1 



^ " rnnpq ^'^rs ^ rpq - 



7^rs ^ rpq ^smn 



diQmdiQndjQpOjQ q I Q^rs ^rmn'^spq Q^mnpq 



(4.10) 
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See Appendix ^ for detailed derivations. 

After a straightforward calculation and changing into comoving time 77 defined hy dt = a drj, we get the leading-order 



4th-order (interaction picture) interaction Hamiltonian in terms of and Qs, for the model (3.5): 



with 



(4.11) 

(4.12) 
(4.13) 



ecQ',Q'si^^Qa^^Q,) 



Q'^ {dQsf + Qsa Q's [OQaf 
+ {OQaf (dQsf + {^^Qa^^Qsf , (4.14) 

where the various coefficients can be found in Appendix ^ It can be seen directly from (4.11 )-( j4.14 ) that there are three 
types of four-point interaction vertices, one involving four temporal derivatives, one involving four spatial derivatives and 
one involving two temporal and two spatial derivatives. This is similar to the case in general single field inflation. Actually 
as we will see below, the four-point functions can be grouped into three types which correspond to three fundamental 
momentum-dependent shape functions. Moreover, there are 4Qcr, 4Qs and IQ^lQs interactions. Thus the corresponding 
non-vanishing four-point correlation functions are {QaQaQaQa), {QaQaQsQs) and {QsQsQsQs)- 



5. The Trispectra 

In this work, we focus on the tree-level four-point correlation functions from direct four-point interactions. In the cosmo- 
logical context, correlation functions are conveniently calculated by using the so-called "in-in" formalism (see Appendix 
^ for a brief review). For our purpose, the tree-level four-point functions are evaluated in the following form 



drj' (0 10(77,) il4(^')|0) 



(5.1) 



where O denotes product of four fields, e.g. QaQaQaQa etc., and H^^ is given in ( paT| ). Although we do not write down 
them explicitly, we should keep in mind that, all quantities in (5.1) are "interaction-picture" quantities, and thus |0) is the 
free vacuum. 

5.1 Four-point functions of the inflaton fields 

Since there is no tree-level two-point cross correlation between and Qg, i.e. (QaQs) = 0, the tree-level four-point 



functions are directly related to the corresponding four-point interaction vertices. The various coefficients in (4.12)-(4.14) 



act as "effective couplings" of four-point interactions. From Appendix^ they are combinations of H, e, c^, Cg, A and 11, 
thus in this work, we treat them as approximately constant. 

5.1.1 (Q.Q.g^Q,) 

There are three types of four-point adiabatic mode self-interaction vertices, as shown in fig.|l| And it is easy to read the 
corresponding contributions according to the Feynman-diagram-like representations in fig jl|: 



= (271)353 (fci234)r. 



24r„i 



-00 077 0^7 0^7 0^7 



ki for short 



where we define K = ki + k2 + ^3 

{Q (tQ aQ aQ a} 2 



US 



(5.2) 



d d 

49^ (-fc3-fe4)7 / ^77' —Gfci (77,, 77)— 0^2(77,, 77)0^3 (77,, 7;) (77,, 77) + Sperms 



{~k3 ■ ki) 



"-1 "-2 



1 + 12- 



.hki 



K 



Sperms 



(5.3) 
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Figure 1: Diagrammatic representations of the four-point adiabatic mode self-interaction vertices. A "red 
dot" denotes derivative with respect to comoving time rj. A blue dot denotes derivative with respect to 
spatial coordinates, or in Fourier space the spatial momentum. A line between two blue dots denotes "scalar- 
product". 



where "5 perms" denotes other 5 possibilities of choosing two out of the four momenta fci • • • fc4 as scalar-products. 



8 (fci • fc2) (fca • ^4) ^ai / d?y'Gfci(r/,,77)Gfc2(?7,,ry)Gfc3(?7,,77)Gfcj77*,r/) + 2perms 



-2I£_ 



(fcl ■ fc2) (fc3 ■ fc4) (^^^^Ihh^Ih<^^^Ilh fj^l 



K 



if3 



+ 2perms 
(5.4)' 

where the other 2 permutations are combinations correspond to scalar-products (fci • ^3) {k2 ■ fc4) and (fci • ^4) (fc2 • ^3)- 
It will prove convenient to define three fundamental shape functions'^: 



Ai{ki,k2;k3,k4, 



K5 ' 



A2{ki,k2;k3,k4) = (-fea • k^) 
A3(ki,k2;k3,ki) = — 



k^kA k^ + kA 
1 + 12-^ +3 



if2 



K 



(5.5) 



It is useful to note two properties of these shape functions: 



• permutation symmetries: 

Ai{ki, ^2; ^3, ^4) is completely symmetric with respect to the four momentum ki, k2, k^ and ^4. A2{ki, fc2; fea, ^4) 
is symmetric under permutations fci ^ k2 and k^ <-> fc4. Similarly, A3{ki, k2; k^, k^) is symmetric under permu- 
tations of fcl fc2, ^3 ^ ^4 and (fcl, fc2) ^ (fc3, fc4). 



• scaling properties: 

Ai (Afci, Afc2; Afca, Afc4) = A^Ai (fci, ^2, ^3, ^4) , 

A2 (Afci,Afc2; Afc3,Afc4) = A^A2 (fci,fc2,fc3,fc4) , (5.6) 
A3 (Afci, Afc2; Afcs, Afc4) = A^Aa (fci, fc2, fc3, ^4) • 

In other words, Ai's have momentum dimension as k^: [Ai] = [k^]. 



'^These three shape functions have been found by several authors in investigation of trispectrum in general single-field models f^,^ |62|,|3l||. Here in 
this work, for later convenience, we do not include the permutations in the definition of Ai's themselves directly, in order to generalize them to the case 
of multi-field case with Ca ^ Cc. The reason is that in single-field case, the shape functions have apparent permutation symmetries, which ai'e broken in 
multi-field models. 
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In terms of these three shape functions, we can write 



14 



{QaQaQaQa)^ = (27r)^5^ (fci234)r. \ ^^[[^ J ' ; ^3 , ^4) 

= (27r)"53(fei234)jY% (-72r,)^i ic.h,c,k2;c,k3,c,h) 



(5.7) 



where we have used the scaling property Ai (cafci, Cafc2, c-dk^, c^k^) = c^Ai {ki, k2, k^, k^). It will prove convenient to 
use the "sound speed dependent" (in the later we denote "cs-dependent" for short) shape functions Ai (cafci , CaA;2 , Cafcs , Cafc4), 
especially when the adiabatic mode and entropy mode have different sound speeds: Ca ^ Ce- Using the "cg-dependent 
shape functions" not only makes the expressions simpler and more symmetric but also make the physical picture clear 
Similarly, we have 



7?4 2. 

{QaQaQaQ<j)2 (Stt) "(5^ (fci234) 6^7 — [^2 (Cafcl , Cafc2 ; Cafes , Cafe4) + 5 pCrms] 

1 U ^a 

{QaQaQaQa)^ = {2n)^^ {k^sd (-2f^<,) — [^3 (cafei , CaA;2 ; €3^3 , Cafc4) + 2 pcrms] 

1 l i ^a 

The whole contribution from (QaQaQaQa) is given by 



(5.8) 



(QaQaQaQa) = ^ (QaQaQaQa), 



(5.9) 



5.1.2 {QaQaQsQs) 

There are six types of adiabatic/entropy four-point cross-interaction vertices, involving two adiabatic modes and two 
entropy modes, as shown in fig.|[ 






.■■k3 



V, 

k4 ' 



V, 

k4 \ 



Qa 



Qa 



ski 




ski 




/ks 



k4 \ 



if 

/ks 



k4 \ 



Qs 



Figure 2: Diagrammatic representations of the four-point adiabatic/entropy modes cross-interaction ver- 
tices. A "red dot" denotes derivative wrt comoving time -q. A blue dot donotes derivative wrt spatial coordi- 
nates, or in Fourier space the spatial momentum. A line between two blue dots denotes "scalar-product". 



= {27rf6\k,23A)r^ 



/^* d d d d 

oo d?? d77 d77 d77 

8\ J2^2 1,2 



i2ff8\ Casern fcr 



3 



(5.10) 



'''in the last line of ( |5.7| ), we abstract a factor rather than in previous works, in order to make the expression more symmetric. The reason will 
become clearer in the following calculations for {QaQaQ sQa) ■ 
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where K = (fci + fc2) + Ce (fca + k^). In terms of the three fundamental shape functions (5.5), (5.1C) can be written in 
a rather compact and convenient form: 

{QaQaQsQ.), = (27r)"(53(fci234)^^(-12r,)Ai(cA,Cafc2;Cefc3,Cefc4) , (5.11) 



where Va and Vs are the dimensionless power spectrum defined in ( |3.23[ ). Similarly, we have (readers who are interested 
in their explicit expressions may refer to Appendix ^ 

(QcrQaQsQs)2 = (27r)"j'^(fci234) J^* ,T — f^^2 (Ca^l , Cafc2 ; Cefca , Cefc4) , 
{QaQaQsQs)^ = (27r)"(53(fc^23^)^i^ {c,h , C,h; C.k^ , C,k2) , 

{QaQaQsQs)i = (27r)"(5^(fci234) J^* ,3* — [A2 (CaA:i,Cefc3;Cafe2,Cefc4) +^2 (Ca^l , CeA:4 ; Cafc2 , Cefes) 

+ ^2 (Cafc2,CeA:3;Cafcl,Cefc4) +^2 (Cafe , CeA:4; Cafci , Cefcs)] , 
(Q.Q.QsQ.)5 = (27r)"(53(fci234)^^^ ^3 (Cafcl, Cafes; Cefca, Cefc4) , 

{QaQcrQsQs)e = (27r)"(5^(fei234) j^* ^3* 1^3 (Cafel , Cefes; Cafc2, Cefe4) +^3 (Cafel , Cefe4 ; Cafcs , Cefes )] • 

(5.12) 

For clarity and later convenience, we group these six contributions into three types, arising from Ai, A2 and A3 
respectively. First we have, 

(QaQ.QsQs)A, = {QaQaQsQs)i = (2^) " <5' (fcl234) (-12r,) Ai (Cafci, Cafczl C,h, C,ki) , (5.13) 

While 



s/4 

= (27r)"(5^(fci234) T^* , r i — f^^2 (Cafcl,Cafc2;Cefe3,Cefc4) H ^^2 (Cefc3, Cefc4; Cafcl , Cafes) 

Qc 

+ T [A2 (Cafcl, CeA:3; Cafes, Cefe4) + A2 (Cafci , Cefc4 ; Cafc2 , Cefes ) 

+A2 (Cafc2,Cefc3;Cafel,Cefe4) + A2 (Cafc2 , Ce A;4; Cafel , Cefe3 )] | , 

(5.14) 

and 

{QaQaQ sQ s) = {QaQaQ sQ s) + {QaQaQsQs)Q 

= (27r)"(53(fc^234)^^^ I ^^f^r^ ^3 (Cafcl, Cafc2; Cefe3, Cefe4) (5 jg^ 

+ 2c2|;2^ t"^^ (Cafcl,Cefc3;Cafc2,Cefe4) + A3 (Cafel , Cefe4; Cafc2, Cefc3)] 

The whole contribution from (QaQaQsQs) is given by 

(Qa{kl,f]*)Qa{k2,V*)Qs{k3,V*)Qs{f^4^V*)) = {QaQaQsQs) + {QaQaQsQs) + {QaQaQsQs) A3 ■ (5-16) 

5.1.3 (QsQ.QsQ.) 

There are also three types of Qs self-interaction four-point vertices, as shown in figj|. {QsQsQsQs) are easily obtained 



by simply changing the corresponding "effective couplings" and Va ^ Vs and Ca ^ Cg in ( |5.7[ )-( |5.8[ ) 
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/ \ / \ / \ 

G. G. G. G. 

Figure 3: Diagrammatic representations of tiie four-point entropy mode self-interaction vertices. A "red 
dot" denotes derivative wrt comoving time r\. A blut dot donotes derivative wrt spatial coordinates, or in 
fourier space the spatial momentum. A line between two blue dots denotes "scalar-product". 



-p4 

{QsQsQsQs)i = (2^)"^^(fei234)— f^(-72r,)Ai (Cefci,Cefc2,Cefc3,Cefc4) 



7?4 Y 

{QsQsQsQs)2 = (27r)"(53(fci234)— f^e^^ [yl2(cefci,ce/c2;cefc3,cefe4) +5perms] (5.17) 

Hi H 



VI , „^ . 1 



{QsQsQsQs)^ = (27r)"(5^(fci234) ^/*^ (-20^) ^ [A3(Cefcl,Cefc2;Cefc3,Cefc4) + 2 pCrms] 



The whole contribution from (QsQsQsQs) is given by 

3 

{Qs{ki,Tj*)Qs{k2,v*)Qs{k3,ih)Qs{k4,V*)) {QsQsQsQs)i ■ (5.18) 



5.2 Four-point Function of the Curvature Perturbation (7?.^) 

As we have addressed before, the scalar-field perturbations Qc and Qs themselves are not directly observable. What we 
are eventually interested in is the curvature perturbation TZ. As has been investigated in [ |57[ ^ |62[ |, this can be achieved 
by writing 



= NaQa* +MsQs* ■ 



(5.19) 



Note that Ms = T-izsMa is in general time-dependent. The four-point correlation function for comoving curvature 
perturbation TZ is given by 



(7^(fcl)7^(fe3)7^(fe2)7^(fe4)) 

-N^ {Qa{ki)Qaik2)Qa{k3)Q,iki))^ +NlN^ [(Q,(fci)Q,(fc2)g.(fc3)Q.(fc4)), + 5 perms] 
+ Nt {Qs{ki)Qs{k2)Qs{k3)Qs{ki)), , 



(5.20) 



where "5 perms" denotes other 5 possibilities of choosing two momenta for Q^, and two momenta for Qs out of the 
four momenta'''. In deriving (5.20), we have used the assumption that there is no cross-correlation between adiabatic 
and entropy modes around horizon-crossing, i.e. (QcrQs), = 0. This is indeed the case if we use the second-order 



perturbation action (3.8) as our starting point of canonical quantization. 



It is also convenient to group the whole contributions to (7?.'') (5.20) into three types, which correspond to three 
different physical origins or more precisely three different types of four-point interaction vertices. From (|5.7|), (5.13 ) and 



i^That is, {Qc{ki)Qa(k-i)Qs{k2)Qs(h4,)), {Qa(ki)Qa(k4)Qs{k2)Qs{ki)), 

{Qa{k2)Qa{k4)Qs{kl)Qs{k-i)) and {Qa{k3)Qa{k4)Qs{ki)Qs{k2)). 



{Q<y{k2)Qa(ks)Qs{kl)Qs{k4)), 
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(5.17), we have 



(7^*)^ = (2^)",53(fei234)j^ WaVt. (-12r,) Ai {c,h,c,k2;c,k3,c,k4) 

+ U^U!vj,Vl (-12r,) Ai {c,kuc,k2; c.k^, ck^) 
^^s'Ps* (-12r^) Ai (cefci, Cefc2, Cefcs, c^ki) + 5 pems} 

= (2^)"j^(fci234) -^^g'^ jr^Ai {c,ki,C,k2;C.,k3,C,ki) 



TnS— Tc Ai (Cafci, Cafe; Ce^a, Cefc4) 



(5.21) 



) Ai (Cefci, Cefe, Cefe, Cefe) + 5 pemS 



where we have used A/^ = T-jisMa and the fact that VsjVa = Ca/cg. "5 perms" denotes other 5 possibilities of choosing 
two out from the four momenta fci, • • • , A:4 as the first two arguments of and the other two momenta as the last two 
arguments of Ai(obviously the differences only arise in A\ {c^ki, Csk2] c^k^, Cefe)). Similarly, we have 



(n^)^ = (271)1^3 (fci234)^^^^ (%A2(Cafcl,Cafc2;Cafc3,Cafc4) 



ns- 



—^-^2 (Cafcl,Cafc2;Cefc3,Cefe4) + A2 (Cefe, Cefe; Cafcl , Cafc2) 

'^e ^a 

[A2 (Cafci,CeA:3;Cafc2,Cefc4) +^2 (Cafcl , Cgfe ; Cafc2 , Cefc3 ) 



-A2 (Cafe, Cefe; Cafci, Cefc4) + A2 {c^k2, C^k^, C^ki, Cgk^ 



T-RS— I ^^2(cefci, Cefe; Cefc3, Cefc4) + 5 perms \ , 



(5.22) 



and 



(7^^)3 = (2.)"^^(fe,234)^ 



-A3(Cafci, Cafc2; €3^3, Cafc4) 



7^5- 



^ f 2"'* ^3 (CaA;i,CaA;2;CeA;3,Cefc4) 



(-»c) 

2c?c? 



[A3 (Cafcl,Cefc3;Cafc2,Cefe4) + Aj, {cjlx, cjft^; cjt,2, 



K<S — 



A3{ceki,Cek2;Cek3,Ceki) +5 pcrms 



(5.23) 



It is convenient to define three new shape functions, which can be viewed as the generalizations of the three funda- 



mental shapes (5.5) 



Al (fci,A:2,fc3,fc3) = Ai (Ca/Cl, Cafe; Cafe, Cafe) + T, 



■ Al (Cafcl, Cafe; Cefe, Cefe 



Ttzs- 



(5.24) 



■ Al (cefe, Cefe, Cefe, Cefe) + 5 perms , 



-17- 



^2(fel, fc2; fcs, ^4) = ^A2(Cafcl,Cafc2; Cafes, Cafc4) 



1 



%^A2 (Cafci, Cafc2; Cefcs, Cefe4) + ^|^^2 (Cefcs, Cefc4; Cafcl, Cafc2) 



+ T [^2 (Cafcl,Cefc3;Cafe2,Cefe4) +^2 (Ca^l , CeA:4 ; Cafc2 , Cefcs ) 



-^2 (Cafc2,CeA;3; Cafcl, Cefc4) +^2 (Cafc2 , Cefc4 ; Cafcl , Cefcs) 



, ^ A2(cefci,Cefc2;Cefc3,Cefc4) +5 perms, 



(5.25) 



and 



-43(fcl,fc2;fc3,fc4) = ^A3(Cafel,Cafe2;Cafe3,Cafc4) 



+ \Tns-] ^ 



2c2c2 
4 



CaCe 



A3 (Cafel, Cafe2; Cgfes, Cefe4) 



+ 773^ [^3 (Cafel,Cefc3;Cafc2,Cefc4) + A3 {c^ki, C^ki, C^k2, C^ks) 



(5.26) 



+ T-JZS— --A3(Cefel, Cefc2; Cefc3, Cefc4) + 5 pcms . 



Note that Ai also depend on parameters Ca, Ce, A, 11 and T-jis- In multi-field models, these parameters enter the definition 
of the shapes. This is essentially different from that in single-field models, we can always abstract shape functions, 
which are functions of momenta only. In terms of these three "generalized shape functions" (5.24)-(5.26), the four-point 
correlation function for comoving curvature perturbation TZ can be recast into a rather convenient form: 



(7^(fcl)7^(fc2)7^(fc3)7^(fc4)) = (27r)"53(fci234)^w^ [-i2t,Ai + e,A2 - n.A^] 



(5.27) 



5.2.1 Trispectrumof 7^ 

In practice, it is convenient to define a so-called trispectrum for TZ: 



(7^(fcl)7^(fc2)7^(fc3)7^(fc4)) = (27r)^5^(fci234)T7j(fci,fc2,fc3,fe 



(5.28) 



From (5.27) we have 



Tn{ki,k2,k'i,ki) 



1 

Ipf 



{-l2T„Ai + QaA2 - n,Az 



2H\ 



(5.29) 



i-UT^Ai + eaA2 - n^As) 



where we have used A/^ = l/2e and Va* = ^ (5^) . At this point, it is useful to note the dimensions of various 
quantities: 

. [r.] = [9.] = [n,] = [i/H% 

. [Ai] = [A2] = [As] = [k% 

• thus, the trispectrum has dimension as [T^] = [H^/k^]. (Recall that the power spectrum has dimension [Pn] = 
[H'^/k^], and the bispectrum has dimension [B-jz] = [H^/k^].) 



To end this section, we would like to make several comments. 
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• In this work we investigate multi-field models with Lagrangian of the form (3.5), which are general function of 
two independent contractions X = Xj and Y ~ XjXj . One may expect at first that the parameter space for the 
trispectra will become more complicated. However, as we have shown explicitly in this work, the final (leading- 
order) trispectrum for TZ is controlled by six "parameters": e, Ca, Ce, A, 11 and Tjzs- Especially, the non-linear 
structure of P{X, Y, (p^) affect the final trispectrum through c^, Cg, A and 11, which are specific combinations of 



derivative of P with respect to X and Y (see (3.7) and (3.10)). On the other hand, this fact would make it more 
difficult to determine the structure of a multi-field model from observations, since the functional forms of P are 
highly degenerate with respect to the parameter space of trispectra'*. 

In the limit of small speed sounds (ca <C 1, Ce ^ 1), at the leading-order, there are three types of "contact four- 
point interaction vertices" according to the "derivatives": vertices with four temporal derivatives, vertices with 
two temporal derivatives and two spatial derivative, vertices with four spatial derivatives (as illustrated in e.g. 
Fig|l|). This is exactly the same as in general single-field models, where three types of contact interaction vertices 
correspond to three different "c^-independent" shape functions. However, the subtlety is that, although in multi- 
field case the interaction vertices can still be grouped into these three types as in single-field case, since there are 
more fields and more speeds of sound, parameters arising from the non-linear structure of the theory enter into the 



definition of shape functions, as in (5.24), (5.25) and (5.26). Thus, in general there is no hope to abstract pure 



momentum-dependent shape functions without involving e, Ca, Ce, A, 11 and T-jis, and thus the so-called shape 
functions in general take the form 

A{ki,k2, fca, ^4; e, Ca, Ce, A, n, T-jzs) . 
This fact makes the analysis of the shapes of non-Gaussianities in multi-field models much more complicated. 



6. Characterizing the Trispectrum 

For the power spectrum, it is easy to plot a curve for P{k) which is a function of a single variable k. However, for higher- 
order correlation functions, or their fourier transformations bispectrum, trispectrum, etc, it is not easy to abstract simple 
quantities to characterize them. While from the point of view of comparing theories with observations, this is obviously 
one of the most important aspects in higher-order statistics of cosmological perturbations. 

In this section, first we analyze the "parameter space" for the trispectrum or general higher-order correlation functions. 
Base on the discuss on "inequivalent momentum configurations", we plot various "shape" diagrams which characterize 
the trispectrum on a 2-dimensional subspace of its whole parametric space. These "shape" diagrams are efficient for 
characterizing the trispectrum (or higher-order correlation functions) and are also convenient for visualization. In the end 
of this section, we map the trispectrum to real numbers, which can be seen as the analogue of non-linear parameters tnl 
and (7nl in the literatures. 

6.1 Parameter Space for the Trispectrum 

In general, the four-point correlation function, or its Fourier transformation, the trispectrum, is a function of four 3D 
momenta: T (fci, fc2, ^3, ^4)- Let us first analyze the "dimension" of the parameter space of T (fci, fc2, ^3, ^4), i-C-, the 
number of independent degrees of freedom in order to identify a "inequivalent momentum configuration" (in other words, 
the dimension of "equivalent class" of momentum configurations). At first, we need 12 real numbers to specify a set of 
four 3D momenta {fci, fc2, fc3, fc4}. However, the 3D momentum conservation (fci + fc2 + fes + fe4 — 0) eliminate 3 
of them, and since the cosmological perturbation is assumed to be statistically homogeneous and isotropic'^, we have 3 
dimensional rotation symmetry SO(3) which eliminates another 3 of them (that is, we call two momentum configurations 
are "equivalent" if they can be related by SO(3) transformation'^). Thus, the dimension of parameter space for trispectrum 
is 6. This argument can be generalized to arbitrary higher-order correlation functions. In general we need 

3n-6, (n>3) (6.1) 

'*This is very different from single-field case, where the large trispectra are expected determined by three parameters Cs, A and 11 (or their functions) 
which are combinations of three derivatives P^xx< P,xxx 2nd P^xxxx- Thus in principle one may strict the functional form P{X) fairly when 
Cs, \ and n have been known. 

"See for a recent investigation on models which break statistical isotropy, where the quantum fluctuations are generally statistical anisotropy, 
and the power spectrum not only depends on fc = |fc|, but also depends on the orientation of k. 

'^Of course there have some discrete symmetry, e.g. permutation symmetry, but this will now affect the dimension of continuous parameter space. 
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Figure 4: Pictorial representation of sets 
of six parameters {^1,^2,^3,^4,^12,^23} or 
{fci, fc2, fci2, 6*1, ^2, ^3} to specify an inequiva- 
lent 3D momentum configuration for trispectrum. 



Figure 5: Pictorial representation of a sef of five 
parameters {ki, fc2, ^3, ^4, fci2} to specify an in- 
equivalent planar momentum configuration for 
trispectrum. 



real parameters to specify a momentum configuration for a n-point function. While case for n = 2 is trivial, it is well- 
known that the power spectrum is a function of one parameter k = \k\. For example, the three-point function or its fourier 
transformation bispectrum takes the form B{ki, k2, k^), which depends only on the modulus of three momenta, which 
are the nature choices of parameters for the bispectrum. In the case of trispectrum, things become much more involved, 
since we need 6 parameters. 

• One possible choice is to choose a set of 6 parameters {ki, fc2, fca, fc4, ki2, fc23} with fci2 = |fei + fc2| and k23 = 
1^2 + ^31- A pictorial presentation is depicted in fig. ^. At this point, it is useful to write down other four 

kij = \ki + kj\ in terms of {/ci, fc2, fca, A;4, fci2, A:23} for later convenience'^ 



fcia = A:24 = y fc^ + fc| + /c| + fc| - fc^2 - ^23 : 2) 
hi = k23 , hi = fci2 • 

Obviously, these six momenta cannot take all real values simultaneously. Actually they should satisfy the so-called 



"triangle inequalities" (remember we have (6.2)): 



h + kj > hj , h + hj > kj , — 1,2,3,4: and i ^ j) ■ (6.3) 

Another possible parameterization is to choose three momenta and three angles: {fci, k2, /c23, ^1, &2, O3} (see also 
fig.H). Similarly, the three angles 61,62, 63 should satisfy: 



0^ + e,>ek, {i,j,k = 1,2,3 and i^j^k), 

h + e2 + e3<2TT. 



(6.4) 



There can be simplifications from the point of view of comparing theory to observations. Since the Cosmic Mi- 
crowave Background is essentially a two-dimensional statistical field, there is actually one constraint among the above 
six momentum parameters. Most significantly, the "planar momentum configuration" is of special importance, since 
in practice one average over numbers of small regions in the CMB, which are approximately planar Thus, in the two- 
dimensional case, we need (2 dof from momentum conservation and 1 from SO(2) rotation symmetry) 

2n - 3 (6.5) 

parameters to characterize a planar momentum configuration for n-point function. For example, for three-point function 
(or bispectrum), the momentum configuration is always a triangle, thus there are 3 independent parameters to determine 



"Remember that since momentum conservation is: fci + fc2 + ^3 + fc4 = 0, it immediately follows A:i3 = |fci + fc3| = |fc2 + fc4| = ^24- 
^"Future LSS (Large-Scale Structure) experiments will bring us valuable information on 3D-configuration of the trispectrum. 
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a configuration of bispectrum. While we need 5 independent parameters to specify a planar momentum configuration for 
trispectrum. How to choose these five parameters properly is somewhat tricky. In this work, we choose the five parameters 
as {fci, k2, fca, /c4, ki2}, which is depicted in fig^. Moreover, these four momenta must satisfy a "planar condition", this 
condition can be expressed as writing k23 in terms of other 5 momenta: 



^23 ^ ykf + kl- 2fcifc4 cos(a + /?) , 

l2 I l2 _ 1.2 i„2 ,1.2 _ 1.2 (6.6) 

with cosa= ^, cos/3= ^. 

2fc4Ki2 2fclKi2 

6.2 Shape of the Trispectrum 

In this section we plot the shape diagrams of the trispectrum, i.e. two dimensional surfaces which capture parts of the 
properties of the trispectrum^^ In this work we plot the shape diagrams with one 3D momentum configuration and one 
2D momentum configuration. 

• For 3D momentum configurations, we would like to choose ki = k2 = ki2 — k23- Thus we are left with two free 
parameters which one may choose as k^ and ki. Actually in this case, it is more convenient to choose two angles: 
9i and 02, which satisfy the inequahties ^ < 9i + O2 < ^ and \0i — 021 < f • Note that we have 



— y k2 + fcfs — 2/02^23 cos 02 , ki — ykf + fcfg — 2fcifc23 cos 0i . 

The shape diagrams of Ai, A2 and A3 are depicted in fig.|[ 

• For planar momentum configurations (see Fig.||), we consider fcs = ^4 = ki2, thus the parameter space becomes 
{ki, k2}. The shapes are depicted in fig.0. 

6.3 Non-linear Parameters 

In practice it is also convenient to define some so-called non-linear parameters, which characterize the typical or overall 
amplitudes of non-Gaussianities^^. For bispectrum which is defined as 

(7^(fcl)7^(fe2)7^(fc3)) = i27rfS''{ki23)Biki, k2,k3) , (6.7) 

the non-linear parameter /nl is usually defined as 

B{ki,k2,k3) = \hdkl,k2,k3)[P{ki)P{k2) + P{k2)P{k3) + P(k3)P{ki)] . (6.8) 

5 



Note that in general /nl(^i, ^2, ^s) is a function of three momenta. ( p.8[ ) is motivated from the "local"-type bispectrum, 
which arises from interaction of the form as local product in real space: 

7^ = 7^g + ^/iJT' {nl - {n,f) , (6.9) 

where 7?.g is Gaussian. In this case, /nl is a pure real number, and the bispectrum is maximized in the limit of one of 



the three momenta going to zero. Although interaction of the form ( |6.9| ) is simple in real space, it is difficult to generate 
large local bispectrum in realistic models. Actually this form of non-Gaussianities are mostly produced by nonlinear 
gravitational evolution subsequent to the horizon-crossing [B2l 93] . 



For trispectrum as defined in (5.28 ), the case is much more complicated. One possible parameterization is 



54 

T{k,,k2,k3, ki) = tnl [P{k^3)P{k3)P{ki) + 11 perms] + —g^^ [P{k^)P{k2)P{k3) + 3 perms] , (6.10) 

25 



^' Since the momentum configuration space for the trispectrum is high-dimensional (6 for 3D configuration, 5 for 2D configuration), the proper 
visualization of the trispectrum is still a challenge. The traditional "shape diagrams" are two-dimensional surfaces which captures parts of the properties 
of the trispectrum, i.e. the projection onto a two-dimensional subspace. 

-^Mathematically, this is essentially to map T(fci, fc2, fca, fc4) to a real number. The "inequivalent momenta configuration" space M for the 
trispectrum T{M) is 6-dimensional. Traditional definitions of non-linear parameters con'espond to choosing one specific point p £ JV[ (i.e. one 
specific momenta configuration), and defining the non-linear paramter(s) proportional to T|pgjVi S M. 
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3 "^0 

Figure 6: Shapes of ^i, A2 and ^3 with momentum configuration ki — k2 — fci2 = ^23 as function of 9i and 62. In this group of 
diagrams, parameters are chosen asca = Cc = 0.1,e = 0.01,A = n = 100 and T-jzs = 1- 



where kij = \ki + kj\. This parameterization of trispectrum is also motivated by the "local-product- type" interaction in 
real space: 



n = n„ 



local \ 
'NL ) 



(6.11) 



2 V iNi. / V g ^ ^, J 25- 
Actually, ( pO| ) itself has no prior convenience for the parameterization of trispectrum in general context. This is justified 



by the fact that, in models with non-canonical kinetic terms, especially with small speed(s) of sound, the leading-order 
non-Gaussianities are originated from the "derivative-coupled"-type interactions, which is obviously not real space "local- 



product"-type interactions described as (6.9) and (5.11). 



Actually, there is a general definition of non-linear parameters for the trispectrum. From (5.29) and the discussion in 



sec. 5.2. 1 , the trispectrum has dimension -p,-, thus one may conveniently define a dimensionless non-linear parameter tnl 
by 



/2 \6'p3 rNL(fcl,fc2,fc3,fc4) 

^ ' '"^ F{k,,k2,k3,k,) 



(6.12) 



where V-jz is the power spectrum of curvature perturbation on large scales which is given in ( [3.29| ) and F{ki, fc2, fcs, ^4) 
is some function of momenta which has momentum dimension k^. Remember that in general the so-called non-linear 
parameters tnl are always functions of momenta. 

As addressed before, the non-linear parameter is essentially a characterization of the typical or overall amplitude of 
non-Gaussianities. Thus one may freely choose convenient F{ki, ^2, ^3, ^4) and momenta configurations to abstract real 
numbers for various purpose. In this work, in order to get a glance of the amplitude of non-Gaussianities from trispectrum, 
we abstract one real number tJ,* from the trispectrum, which are defined respectively 



23. 



(6.13) 



^Tjlj* has also been used in a recent investigation of the trispectrum in general single field models | |3l| ]. 
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Figure 7: Shapes of Ai, A2 and A3 with 2D momentum configuration ki — k2 = k\2 as function of k\ and ^2. In this group of 
diagrams, parameters are chosen as Ca = Cc = 0.1, e = 0.01, A = 11 = 100 and Tizs = 1. 



in which "regular tetrahedron" denotes the "regular tetrahedron" momenta configuration for trispectrum: ki — k2 = 
k3 = = ki2 — k23 = k (see fig.H). Here we conveniently choose F = k^. One can of course define other convenient 



non-linear parameters by choosing specific momentum configurations and F in (6.12) for particular purpose. 



As has been discussed in the end of sect. 5.2. 1 , in multi-field case, the shape functions intrinsically involve parameters 



Uke Ca etc. Thus, even the momentum configuration is fixed, the trispectrum or the corresponding non-linear parameters 



defined as in (5.13) is complicated. In this work we give the expression in the limit Ca <C 1 and Ce ^ 1. We have 



rth 
'NL 



rp2 '^a 

^ns — 



(ti +T2+ T3) 



(6.14) 



where 



T2 



T3 



3c2 (54c2a2 - ij2e(3A + lOH)) 



13 (-H^e 



512i?4g2 



256c? i72e 



rpl 



cl (c: 



F 



ns' 



9ca {H^e + 34 {-Icl + 2c2) A) 



8(Ca + Ce)5i?2e 
2 ' 5CaCe + 1 C^^ {-Hh + 3cjX) 



4cf(Ca 



)5i72, 



515 
8192c? 



TZS' 



9Ca(-13c2 + 4c2 



rTi4 



'^^ 1024c5 

13Ca(c2 - 4) 



Ca (-5c2 + 6c2) (Scf + 25c3ce + 43c2c2 + 25cac3 + 5c^) 

64c4(Ca+Ce)5 



256cJ 
515c? 



ns' 



{cl-cl) 



2048c9 



(6.15) 



7. Conclusion and Discussion 

In this work, we studied the inflationary trispectrum in general multi-field models with scalar field Lagrangian of the form 
P(X^'^, (j)^), which is the generalization of multi-field /c-inflation and multi-field DBI models. In our general framework, 
we expanded the perturbation action up to the fourth-order, and calculated the four-point correlation functions for adiabatic 
and entropy modes, and also the trispectrum of the curvature perturbation on superhorizon scales. 

We have shown that the perturbations for adiabatic and entropy modes are enhanced by 1/ca and l/cg respectively, 
where Ca and Cg are propagation speeds of adiabatic and entropy modes respectively and in general c^. In this work 
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we focus on the trispectrum from four-point interaction vertices. In the two-field case, we have shown that there are 
three non- vanishing four-point correlation functions: {QaQaQaQa), {QaQuQsQs) and (QsQsQsQs)- In the leading- 
order, all four-point correlation functions and also the final trispectrum of the curvature perturbation can be grouped into 



three types, which are deformations and permutations of three fundamental shape functions (5.5). These three types of 



contributions arise from three types of four-point interaction vertices. Thus, our result can be seen as the generalization 
of the previous results in general single-field models ^ 
models 



32 1 and also the very recent investigations in multi-DBI 



In single-field models, one can always abstract momentum-dependent shape functions and put all other parameters 



as overall pre-factors. However, as we have seen in three "generalized shape functions" ( 5.24 )-( 5.26 ), parameters such as 
Ca, Ce, A, n and Tns enter the definitions of these shape functions, and make the momentum dependence of these shape 
functions complicated. In this work, after a general discussion on the parameter-dependence of the trispectrum, we plot 



two set of shape diagrams of ( 5.24 )-( 5.26), including 3D and 2D momentum configuration with fixed parameters as Cg etc. 



However, one expects that there should be better characterization of the trispectrum and the shape functions. We hope to 
get back to these issues in the future. 

In this work we focus on the non-Gaussianities arising from the interactions among quantum fluctuations. While 
interactions not only cause non-Gaussianities, they also cause quantum loop-corrections. By collecting signatures of both 
quantum loop corrections and non-Gaussianities, we will obtain a more sensitive test of the physics of inflation [0. 

In this work we focus on the primordial non-Gaussianities which are evaluated around the sound-horizon crossing. 
While detectable non-Gaussianities can also be produced when the curvature perturbation is generated from the entropy 
perturbation(s) on superhorizon scales, or at the end of inflation, or during the complicated reheating process. Since non- 
canonical kinetic terms can arise naturally in string theory inspired models, cosmic string effects should also be considered 



[114, 1151 



Another issue should be address is that, in this work (and also previous works on non-Gaussianities in general multi- 
field models), the cross-correlation between adiabatic mode and entropy mode around horizon-crossing is assumed van- 
ishing: {QaQs) ^ ~ 0. This is good approximation in slow-roll inflation with canonical kinetic term, which corresponds 
to^~ 



in (3.25) and the background strategy is straight. In this work we also assume ^ 
investigate the case where ^cannot be neglected and (QaQs)^ + [ [To^[lO^ [To^. 



0, while it is interesting to 
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A. A Brief Review of "In-in" Formalism 



A.l Preliminaries 



The "in-in formalism" (also dubbed as "Schwinger-Keldysh formalism", or "Closed-time path formalism") [110 111 



112 1 is a perturbative approach for solving the evolution of expectation values over a finite time interval. It is therefore 
ideally suited not only to backgrounds which do not admit an S-matrix description, such as inflationary backgrounds. 

In the calculation of S-matrix in particle physics, the goal is to determine the amplitude for a state in the far past j?/') 
to become some state \'\\}') in the far future, 

(V''|5|V'> = (^'(+oo)|^(-cx3)). 

Here, conditions are imposed on the fields at both very early and very late times. This can be done because that in 
Minkowski spacetime, states are assumed to be non-interacting at far past and at far future, and thus are usually taken 
to be the free vacuum, i.e., the vacuum of the free Hamiltonian i/o- The free vacuum are assumed to be in "one-to- 
one" correspondence with the true vacuum of the whole interacting theory, as we adiabatically turn on and turn off the 
interactions between t = — oo and i ~ +oo. 

While the physical situation we are considering here is quite different. Instead of specifying the asymptotic conditions 
both in the far past and far future, we develop a given state forward in time from a specified initial time, which can be 
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chosen as the beginning of inflation. In the cosmological context, the initial state is usually chosen as free vacuum, such as 
Bunch-Davis vacuum, since at very early times when perturbation modes are deep inside the Hubble horizon, according 
to the equivalence principle, the interaction-picture fields should have the same form as in Minkowski spacetime. 

A.2 "In" vacuum 

The Hamiltonian can be split into a free part and an interacting part: H = Ho + Hi. The time-evolution operator in the 
interacting picture is well-known 



U{r]2,r]i) = Texp 



dt'Hniv'] 



(A.l) 



where subscript "I" denotes interaction-picture quantities, T is the time-ordering operator. Our present goal is to relate the 
interacting vacuum at arbitrary time to the free vacuum |0i) (e.g., Bunch-Davis vacuum). The trick is standard. 

First we may expand |ili(?7)) in terms of eigenstates of free Hamiltonian Hq, |rii(7/)) = (ni|f7i(7/)), then we 

evolve \^i{ri)) by using (A.l) 



(A.2) 



n>l 



From ( A^), we immediately see that, if we choose 772 = — oo(l — «e), all excited states in ( A^) are suppressed. Thus we 
relate interacting vacuum at = —oo{l — ie) to the free vacuum |0) as 



|a(-c^(l-ze))) = |0i) (Oi|r!i) 
Thus, the interacting vacuum at arbitrary time 77 is given by 

|VAC,in) = \ni{ri)) ^ C/(r;, -cx)(l - ie))|17i(-cx)(l - ie))) 



(A.3) 



Texp —i 



dTj'Hiiirj') |0i) (Oi|r!i 



(A.4) 



A.3 Expectation values in "in-in" formalism 

The expectation value of operator 0{ri) at arbitrary time 77 is evaluated as 
(VAC,in|0(77)|VAC,in) 



(VAC,in|VAC,in) 



= 0, 



T exp i 



D(l+ie) 



drj'Huiv) Oi(^)Tcxp -z 



-00(1- 



dr]'Hii{r]') 



(A.5) 



where T is the anti-time-ordering operator. 
For simplicity, we denote 



— cx3(l — ie) = — cx)^ 



-00(1 + ie) — —00 , 



(A.6) 



since, e.g., —00+ has a positive imaginary part. Now let us focus on the time-order in (A.5). In standard S-matrix 
calculations, operators between (0| and |0) are automatically time-ordered. While in (A.5), from right to left, time starts 
from infinite past, or —00+ precisely, to some arbitrary time ry when the expectation value is evaluated, then back to — cx)~ 
again. This time-contour, which is shown in Fig.||, forms a closed-time path, so "in-in" formalism is sometimes called 
"closed-time path" (CTP) formalism. 



Im 



-oo(l-ie) 



-oo(l+ie) 

Figure 8: Closed-time path in "in-in" formalism. 



Rer] 
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A.4 Perturbation theory 



The starting point of perturbation theory is the free theory two-point correlation functions. In canonical quantization 
procedure, we write a scalar field as 

(/>fe(f?) = w(fc,7?)afc + M*(fc,fy)alfe , (A.7) 

where u{k,ri) is the mode function for (t)k(ri) (in practice, Uk{ri) and u^(?7) are two linear-independent solutions of 
equation of motion for which are Wroskian normalized and satisfy some initial or asymptotic conditions ). 

The free two-point function takes the form 



(A.8) 



with 



In this work, we take (A.8) and (A.9) as the starting point. 
Now Taylor expansion of (A. 5) gives 

• Oth-order 



(0) 

0{v)) = (0,|Oi(77)|Oi 



1 st-order (one interaction vertex) 



(1) 

0{r])) =2 Re 



(A.9) 



(A. 10) 



(A.ll) 



2nd-order (two interaction vertices) 



(2) 



0{t])) =-2 Re 



00 + 



di]' 
dr{ 



dv" Oi 



00 + 



dii" (Oi 



d,{i^)H,,{r^')HAil") 



Oi 



(A. 12) 



Here in this work, since we are considering four-point correlation functions from four-point vertices, ( A.ll ) is needed. 



B. Various Expansion Quantities and Useful Relations 

The exact form 
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2P 



+ 



- U ( yU vKL yMN I yKL yMN yU I yMN yU vKL\ 

:^,{IJ){KL){MN) [^2 ^1 ^1 +^2 ^1 ^1 +-^2 ^1 ^1 ) 



-P 



{IJ){KL)M^ 



-P 



(I.J)KL^ 



, 1 fp ( yI'T vKL yMN yPQ\ ^ A P nP ( yI-J yKL Y^iN\ 

+ -^\P,{IJ){KL){MN){PQ){-^1 Ai Ai ) + ^P,{IJ){KL){MN)PQ (.A^ A^ A^ ) 

+&P,{i.j)iKL)MNQ^'Q''{xi'x^^) + ^P,ii.,)KLMQ''Q''Q'" xi' + PjjklQ'Q'Q'^Q''] , 



P3 — P,{IJ}X3 



IJ 



P,{IJ){KL) Xl^'^X^'"' + P^{IJ)kQ^ X2 



KyU 



I ^ p yU vKL yMN 

+ g-r,{/,7)(A'L)(MJV)^l ^1 ^1 



(B.l) 



(B.2) 



-P 



(ij){kl)mQ'^' X{'' X^^ + -P,(ij)klQ^Q^X(''' + -PjjkQ^Q'Q^ , 
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and 



P2 = P^iij)\-^d,Q'd,Q' + \ 



(q'Q^ - 2 ^[.'iVja.g'^)) - 4ai 4^Q^) + (3a? - 2^2) d^l^i 



1 

+ 2 
+2P 



P 



{I.J)K 



PjjQ'Q 



+ i (3a? - 2a2) (q^Q^ - 2 ^^^a'/Ji^.Q'')) + (-4a? + Baias) ^^^Q-") 
+ \ - 12a?a2 + 3a2) (/jq^o > 

+ (3a? - 2a2) 4>o Q'^'^ + (-2"? + 3aia2) ^i^i , 

Q'Q' - 2 ^'^jNld.Q^A - 4ai 0[,^Q'^) + (3a? - 2a2) ^^^j} 



In this work, at leading-order 



X 



ij 



X. 



IJ 



2 ■^-—d,Q'd,Q-' + -Q'Q 



(B.3) 



(B.4) 



(B.5) 



C. Interaction Hamiltonian 



C.l Derivation of interaction Hamiltonian 

In this appendix, we show the explicit steps to get the interaction Hamiltonian. We start from the following leading-order 
Lagrangian (2nd, 3rd and 4th-order respectively): 

2 2a^ 



2 2a^ 

■^4 — ^mnpq QmQnQpQq ^©mnpg QraQndiQpQiQq ~\~ ^ ^^ninpq diQmdiQndjQpOjQq . 

The conjugate momentum of Q„i is simply 
dC 



(C.l) 



dQm 



3 



1 



(C.2) 



2 — Mtfti, -V -V ^ 2 2^^^^^^^^^^^^^ ^^jnnpqQjiQpQq ^ 2 



^jnnpqQji^iQp^iQq ■ 



In order to get the Hamiltonian, we must solve Q,„ in terms of its conjugate momentum Tr^. In general, it is a complicated 
task. However, this can be done perturbatively. We make the ansatz: 



Qm — -^mn ^' ^mnp '^n'^p + -^mnpq T^nT^pT^q + OijT ) 



(C.3) 



where A*-*' 's are of order unity quantities, which we should determine in the following. Note that in ( |C.3| ), for our purpose, 
higher-order terms are not needed. 



(C.4) 



In order to determine A'-*''s in {23), the strategy is that, plug (C3) into ( |C.2| ), and solve A'-*''s order-by-order 

3„ 
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At first-order, in order to satisfy dC4h, we have 



this gives 

^mn = '^mn • (C-6) 



At the second-order, (C.4) gives 

^\pq ^p^q = ^l^^mnlf^sm^np^lr '^p'^r + ^ "inll^ sm^iQ n^iQ I ■ (C7) 

(2) 

Fortunately, we do not need to solve Xnpq explicitly, the above relation is enough for our purpose. Similarly, at the 
third-order we have 



\(3) 

^spql '^p'^qn 



(C.8) 

9 ■ 



— '^^mnl^ms^np^\rs''^p''^'r''^s mnpq^ms^nl ^pr ^qs ^I'^^r'^s ~^ 2^2^™^P1^™s^nl '^l^iQpdiQ 

The fourth-order part of the Hamiltonian density is 

"^4 = TTm ^mlipq TTnTTpTTg - £(4) , (C.9) 

where £(4) is the fourth-order part of the Lagrangian in terms of Q,„ and 7r„i (note that £2, £3 and £4 all contribute to 
-^(4)): 

£(4) = ICmn ^iip ^mrst '^p'^r'^s'^t + '^^mn '^npg'^mrs '^p'^q'^r'^s + "^^mnl ^inpq ^rir ^Is '^p'^q'^r'^s 

— T , A'^' -k -k Fi d b o, a(^) A^^^A^^U'^^ tt tt tt tt cc lot 



Thus, 



2a2 

^ 2^'nmpq ^rnr^ns '^r'^ s^iQ pOiQ q -\- ^mnpq diQ radiQ ndj Q pOj Q q . 

1 3 

^4 = - 2^mn >^npq>^mrs T^pT^qT^rT^s ~ -j^mnl ^'i^pq ^nl^ls^ TrpTTgTTrTTs 

~t~ ^ 2^ ^^^^ ^rnpq '^p'^q^iQn^iQl ^nuipq ^rnr^ns ^pu ^qv ^r'^s'^u'^v (C-11) 
^" 2^^^^pq ^mr^ns ^r^s^iQp^iQq ^ "^^mrLpq ^iQm^iQn^jQp^jQ q • 

C.2 Interaction Hamiltonian in tlie interaction picture 

Now we go to the "interaction picture", where all fields and momenta are "free". The free (quadratic) part Hamiltonian 
density reads 

^0 = '^In^mli'^n - ^(2) = 2^mn '^\n'^\i + l^^mndiQ^^diQ^ . (C.12) 

Here we explicitly write superscript "I", which denotes interaction picture quantities. In free theory, is related with 
ttJjj by (remember that TiJ, = Ho) 

Qm — "o i~ ~ ^mn^n i i'^- m ~ ^mnQn ■ (C.13) 



Now our task is to plug ttJ^ back into the fourth-order Hamiltonian ( C.9 ), to get its "interaction picture form" Ti^ . After a 
straightforward but rather tedious calculation, we get 



9, 



1 A, A, „ ^. „ ^, /l^ 3 



7^4 QmQnQpQq g^'s ^rmn ^spq ^mnpq + ^2^ niQ n^iQ p^iQ q ^^Qjnnpg /^^^ ^rpq^smn 
H J^^iQm^iQn'^jQp^jQq ( Z^rs '^rmn spq ~ -Z^rrm 



(C.14) 



We should keep in mind that the above Hi is the fourth-order interaction Hamiltonian in the "interaction picture". In this 
work, we work in operator formalism in the interaction picture, thus, we use ( C.14 ) as our staring point. 
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D. Coefficients in Hi 

The various coefficients in the 4th-order interaction Hamihonian density ( |4-.12 ), ( 4.13 ) and ( 4-. 14 ) are 



_ -ij2e(3A + ion) + MX'^cl 

^ (7 — 



ct (7 - 4c2) 2 + 4 (c2 + 4) + cl (-13 - 20c2 + IGc^) 



16i/2ec2d 



1 



4i?4g2c4c2 



[H-e (-1 



=e) ' + 3Ac^ (-7 + 4c2) + H\cl (5 - Sc^ + Acf) 



+4 {iH\ + {-8H^e + 6A) cj + ^Hhcj)] , 



Hhcl+Hh(~2 + a) 



3X4 {-2 + 4) 



SH^eclcl 

-24 + eg (2 + 9cg - 6r4) + 4 (-14 + 15cg - 40^ 



{-1 + 4) (4 + 24 i-i + cl)) 



2me44 



-i±4 



_4 (-2 + 4)' + 4 {-A + 34) 



4i?2g 



(D.l) 



(D.2) 



(D.3) 



E. Explicit Expressions for {QaQaQsQs)i 

{Q crQ aQ sQ 2 
= (2^)3<53(fci234)(-23fi) 



d d 

iQasi I di^ {-ks ■ ki) —GkAv*^v)—^k2{'ri*:V)Fk3{'ri*,'n)FkA'n*^'n) 

— oo 



drj 



drj 



(2^)^J^(fcl234)e,, i-k3-ki) 



H^4 klkl 



1 



12c2fc3fc4 3Ce(fc3 + fc4) 



K 



(E.l) 



(2^)3<53(fci234) (-25R) 



d d ■ 
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K 
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{Q (tQ (tQ sQ 4^ 

= {27Tf6^ku34){~2n) 
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drj 



(2^)^,5^(fcl234)e, 



{-k2 ■ fc4) fc2fc| ( 1 



di] 

l2CaCek2ki 3 {Cak2 + Cek^) 



X2 



K 



3 perms 



(E.3) 



The other permutation is (fc2, ^3), (^i , ^a) and (fci , fc4) combination. 
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{Q aQ aQ sQ 



s/5 



= (27r)^5^(fei234)(-23?) 



= {27rY6%ki234)^o 



(fel • fe2) (fes • ^4) 



K 



(E.4) 



12c^c^ fc- 3 

1 + ' ' + -pCaCe (Ce (fci + fe) ^3^4 + Cafcifcs (fcs + ^4)) 
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